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Abstract 

We consider a system consisting of a planar random walk on a square lattice, sub- 
mitted to stochastic elementary local deformations. Both numerical and theoretical 
results are reported. Depending on the deformation transition rates, and specifically 
on a parameter rj which breaks the symmetry between the left and right orientation, 
the winding distribution of the walk is modified, and the system can be in three 
different phases: folded, stretched and glassy. An explicit mapping is found, leading 
to consider the system as a coupling of two exclusion processes: particles of the first 
one move in a landscape defined by particles of the second one, and vice- versa. This 
can be viewed as an inhomogeneous exclusion process. For all closed or periodic 
initial sample paths, a convenient scaling permits to show a convergence in law (or 
almost surely on a modified probability space) to a continuous curve, the equation 
of which is given by a system of two non linear stochastic differential equations. 
The deterministic part of this system is explicitly analyzed via elliptic functions. In 
a similar way, by using a formal fluid limit approach, the dynamics of the system 
is shown to be equivalent to a system of two coupled Burgers' equations. 

1 Introduction 

Random walks are fundamental objects arising in probability. Also they are of 
primary importance in various fields of physics, especially with regard to polymers 
[jj], H and biology. For instance, planar random walks can be used as a representation 
of DNA coding, since the sequence of the four different kinds of codons (A,G) for 
purines and (T,C) for pyrimidines can be considered as a random walk on a square 
lattice: as a rule, (G,C) code the upward and downward jumps, whereas (A,T) code 
the left and right steps |^, ||. It seems therefore interesting to consider random 
geometrical objects as complex systems, and to submit them to some dynamical 
principles, the goal being to develop methods and tools which hopefully might be 
used to tackle more realistic models. 

In this context we will analyze the evolution of an arbitrary sample path Cn of 
length N, generated by a simple random walk in the square lattice Z 2 , and subject 
to local transformations. This stochastic object has a rich structure, plays an 
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important role in probability theory and lends itself to sufficiently wide but non 
trivial results. 

At time t = 0, Cjv is given, and we assume it has been uniformly generated. This 
means precisely that each successive jump (up, down, left and right) building Cjv 
is selected with the same probability 1/4. Eventually Cjv can be constrained to 
be closed or to have fixed extremities. Once the initial configuration is defined, 
the system evolves according to the four local pattern transformations depicted 
in Figure 



1.1 



Only a single point of the walk can be moved at a time, with the 



constraint that no link be broken (i.e. the walk remains always connected). 
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Fig. 1.1: Pattern transition rates. 



Geometrically, these patterns can be expressed as 

'left bend Ml, right bend M3, 
vertical or horizontal fold M2, 
[straight (-►-►) M4, 
and the following local distortions can occur: 



f 



'Ml -> A/3, with rate A+, 

M3 -►Ml, with rate A", 

rotation of M2 of angle ± \ , with rate 7 ± . 

Hence we have defined a global Markovian continuous time evolution of the system, 
with exponentially distributed jump times, the state space of the underlying Markov 
chain being the set of 4 N sample paths (or curves) Cn introduced above. 

This model is somehow a kind of discrete analogue of the Rouse chain jll| , which is a 
popular model for polymer dynamics. There, each point of the chain is harmonically 
bound to its nearest neighbor, and move randomly in space. Some interesting 
statements can be made concerning the winding properties of such chains in 2d 
Jl2| . In this respect, from a probabilistic point of view, we keep in mind that 
winding variables are for Brownian curves, and they give rise to striking limit laws 
under convenient scalings (see e.g. || f|, |^|). 
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The paper is organized as follows. Section 2 presents the basic numerical and qual- 
itative results, which rely on ad hoc discrete event simulation experiments together 
with a convenient graphical interface. In section 3, we propose several possible ways 
of coding the system with the related probabilistic descriptions (this section can be 
only glanced at without too much damage) . The last part of the paper (section 4) 
contains the main quantitative results: relationships to exclusions models, scaling 
and thermodynamical limit, fluctuation analysis via theorems of central limit type. 



2 Numerical experiments 
2.1 Observations 

We shall begin our study with several basic numerical observations. The model is 
purely stochastic and hence well-suited for Monte-Carlo simulations. These have 
been performed with the help of a graphical interface shown in Figure 2.2, which 
facilitates the exploration of the relevant parameter range, together with the display 
of the main regimes and phases of the system. 
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Fig. 2.2: Platform for numerical experiments. 



Several parameters have to be tuned: the number of steps N; the relative position 
of the extremities of the walk (D x ,D y ); the boundary conditions, which either can 
be defined to let the end points move independently, or can be fixed or to tied 
by some periodic boundary conditions; finally, the time constants associated with 
the elementary transformations, t\ with 7, T2 and T3 with A + and A - . Any walk 
fulfilling these conditions is randomly generated at time t — and then evolves 
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stochastically, with a movement depending on the rates and boundary conditions 
given above. After each new event, time is incremented by an amount inversely 
proportional to the number of all possible moves weighted by their respective rates. 
Interesting things happen when we break the chiral symmetry by imposing a de- 




Fig. 2.3: Picture of a random walk of N = 5000 steps, showing the phases of 
the system for several values of r\. Each colored segment represents 
1000 steps, (a) r\ = (the basic scale), (b) r] = 5 (scale=l). (c) 
77 = 12.5 (scale = 1/6). (d) rj = 250 (scale = 1/2). 



tuning between A + and A , in a proportion of order 

.A+-A- 



V 



N- 



A+ + A- 



(2-1) 



For closed walks, four different situations can roughly be observed (see Figure 2.3). 



( a ) V i5 1- The initial configuration belongs to the equilibrium set of typical 
configurations, only fluctuations are altered by the finite value of r\. 

(b) 1 ^ V ^5 6- The system reaches an equilibrium which still corresponds to a 
random walk, the fractal dimension remaining equal to two, but a macro- 
scopic circular drift is observed, yielding a sort of smoking ring. 



2.2 Brownian windings 
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(c) 6 < rj < 50. The smoking ring gets stretched, and the elementary links 
becomes aligned over long distance. Fractal dimension shrinks to one, with 
the apparition of a long range order. Rotational invariance is broken. 

(d) Tj > 50. The system is not able to reach its equilibrium. This typical con- 
figuration (out of equilibrium) exhibits an intricate hierarchical structure of 
bubbles. Smaller bubbles get evaporated into bigger bubbles. Time constants 
associated to these mechanisms grow exponentially with the size of the bub- 
bles. Therefore, the final configuration corresponding to one bubble is never 
reached in the thermodynamic limit. We will refer to this non-equilibrium 
phase as a glassy phase. 

2.2 Brownian windings 

Some macroscopic random variables of interest can be constructed in order to be 
able to follow numerically the evolution of the system. First of all the total number 
of patterns Ml, M2, M3, M4 is a set a variables which can be used to distinguish 
between a folded and a stretched phase of the system. All these number are expected 
to be fairly distributed in the folded phase although in the stretched phase we expect 
the pattern M4 to be in majority. In order to express mathematically the curling 
of the system, we consider variables related to the winding properties of planar 
Brownian curves, which are defined as follows: with each point in the plane is 
associated its winding angle 0, scanned by the random walker around this point. A 



Fig. 2.4: Winding sectors denned for a closed random walk (a). Color gra- 
dient representing sectors for a walk of N = 10 6 steps: in red and 
blue, respectively, positive and negative winding sectors; in green, 
the null sector (b). 

limit law for this variable has been derived by Spitzer for Brownian curves. Actually, 
assuming the length of the curve is set to 21 — Na? when N — > oo and a — > in 
the Brownian limit, the winding angle 0(£) of an arbitrary point has the asymptotic 
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For a closed walk, the value taken by 9 are limited to 2nir, where n € Z represents 
the winding number associated with the point under consideration. 

For each n, the set of points with the same winding number n form a winding 
sector, whose arithmetic area is a random variable denoted by S„(£). Under the 
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Fig. 2.5: Evolution of the total algebraic area (a,b,c). First transition with 
condensation in the first winding sector (a) N = 1000 steps. Un- 
folding transition with metastable states for N = 2000 steps (b). 
Curves (c) and (d), drawn in a logarithmic scale for N = 10 6 , rep- 
resent respectively the slow dynamics and the distribution of the 
S n s. 

above mentioned Brownian limit for closed curves of length t, we have M 



E[S n (£)] = 




2.2 Brownian windings 
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In a similar way, the total algebraic area enclosed by the Brownian curve is defined 
as 



riGZ 



and its distribution is given by Levy's law | 
lim P (A{£) = 2s£) 



2 cosh (its) 



The variable S n (£) is indeed enough for a complete characterization of the phases 
of the system. 





(a) 



(b) 




(c) 



Fig. 2.6: The limit arithmetic area of the first winding sector when r\ varies, 
rescaled by N (a) and by N 2 (b). The dashed line in (b) is computed 
from part [O^ giving the critical value r] c = 2ir separating the folded 
phases from stretched one. Scaling with iV of time-constants for 
the transient regime (c) for r\ = 10; the fit gives a value 2.01 for the 
exponent. 
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3 Sequence coding and generators: two formal approaches 



2.3 Slow dynamics 

The study of winding sectors and of associated variables is especially well adapted 
to describe the evolution of the system, since the presence of curls has a direct 
impact on the S n distribution. A small but finite value of rj < 1 will show itself 
by the existence of a shift in the distribution of S n together with an unbalancing 
between positive and negative winding sectors. When 1 < ry < 6, the distribu- 
tion of S n condenses into either Si or S-i_. However, the behavior of Si or S-i 
with regard to scaling factor does not change, and these variables still scale like 
N. The system reaches its stationary state after a transient regime characterized 



by a single time constant (Figure 2.5a). For 6 < r\ < 50 the system get stretched, 
an unfolding transition occurs. Si or S-i scale now like TV 2 , and the fact that the 
walk is taut get reflected in the distribution of the motifs Mi, M2, M3, M4. When 
r] > 50 we obtained a glassy phase. This is related to the apparition of a hierar- 
chy distribution of meta-stable configurations. The system evolves slowly to the 
rate of bubbles evaporations, small bubbles collapses and produce bubbles of bigger 
size. As a consequence, the transitory regime is completely different. For small size 
systems (see Figure |2.5fc), it is observed that the total algebraic area increases by 
successive steps. These steps correspond to intermediate metastable states, con- 
sisting of bubbles of increasing size. The associated time constants behave roughly 
exponentially with the size of these bubbles. When N ~> 00 (see Figure pTqc), 
a continuous spectrum of time constants is obtained and the convolution of these 
dynamical effect corresponding to different scale end up in a slow dynamical grows. 
The total algebraic area increases logarithmically with time. We observe also (see 



Figure 2.5d) that the distribution of S n seems to have a limit characterized by the 
absence of negative (or positive) sectors, i.e. strictly zero for negative (or posi- 
tive) index n, together with a scaling exponent around 1. Indeed this sequence of 
distributions seems to behave like n _1 , instead of n~ 2 for 77 = 0. Also the glassy 
transition is clearly of first order, with coexistence of a liquid phase (part of the 
walk which remains folded and disordered) and a solid glassy phase. The param- 
eters corresponding to temperature and magnetization can be defined by analogy 
with standard spin glasses. They can be tuned independently, as we shall see later 
on, by letting vary 77 (which is roughly equivalent to the external magnetic field) 
and 7 (pertaining in some sense to the temperature). 



3 Sequence coding and generators: two formal approaches 
3.1 A group decomposition 

The dynamics of the system is Markovian. Its description involves a configuration 
space {a}, a probability measure P a and the generator G a fj of the evolution op- 
erator. Since the elementary transformations are local, and associated with each 
point of the random chain, G may be expressed as a sum of operators 

N 

Ga/3 = y^ffa/3(^), 
i=l 



3.1 A group decomposition 
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Fig. 2.7: First order glassy transition, with coexistence of a solid phase (on 
the edges) and a liquid phase (at the center). The scale is divided by 
a factor of 2 in (c) and of 4 in (d). The blue gradient indicates that 
all windings are negatives. A complete sequence of blue corresponds 
to an increase of 12 winding numbers. 

Here it is to be understood that each g(i) operates on a real space, endowed with 
the base {u> ai (l) ® uj a2 (2)... <g) uj aN (N)} defined tensorially with respect to sites 
i, where each on stands for possible configurations at position i. Therefore, the 
sequence {ai, ...,on} will generically denote a given configuration of the system. 
Then g has the following representation in the base of local operators 

g(i) = \ + a + (i)(l-a-(i)) + \-a-(i)(l-a + (i)) 
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where a + (i),a (i) , R(i) , R 3 (i) are the elementary operators corresponding to the 
transformations given in the introduction. They enjoy the following properties: 

a ± {if = 0, 

a ± (i)a^{i)a ± {i) = a ± (i), 

a ± (i)R(i) = R(i)a ± (i) = 0, 

R 5 (i) = R(i). 

The fact that the transformations are local is illustrated by the commutation prop- 
erty 

h{i)h(i + 2) = h(i + 2)h(i). 

Many representations can be defined to encode system configurations. For example, 
we could take the sequence 

£! = {zi = x + iy, z 2 , z N }, (z t eZ + it) 

formed by the affixes of the successive points in the plane, with the constraint that 
each point is separated from its neighbors by a single link. Another possible choice 
is 

fii = {v x = z 2 - zi,v 2 = z 3 - z 2 , ...,wjv-i =z N - zjv-i}, {vi e {l,i,-l,-i}), 

the set of successive links, which corresponds to the tangent map to Oo- In this case, 
the center of gravity, i.e. the average over all point positions z* = (z\ + .. + zn)/N, 
defines equivalent classes. Then taking the tangent map to fii yields the sequence 

^2 = {Ul =V 2 - Vi,U 2 , Uat_2 = V N -i - V N ^ 2 }, 

which represents the patterns defined in introduction. In that case, the drift 

v* = (v\ H v N -i)/N 

depicts the separation between extremities, and may serve to define different equiv- 
alent classes. We will see later that, depending on the boundary conditions , z* 
and v* can be decoupled from the dynamics. In particular, this is true for periodic 
boundary conditions. Indeed ^2 should be the natural representation to express 
the dynamics of the system, since elementary transformations are performed on 
the patterns defined at each point of the chain. However, although any possible 
transition at a given point depends solely on the pattern of the point (for example 
a A + -transition at point i can take place only if the walk performs a left bend at 
that point), the result of this transformation involves three successive patterns (for 
the event A + , the pattern at i becomes a right bend, and both adjacent patterns at 
i—1 and i + l are modified). On the other hand, in the Oi representation, although 
each transition is conditioned by a pair of links, the result does modify only this 
pair. 

We choose now Sli to express G properly. As a site i of the chain has four pos- 
sible configurations, depending on the orientation of the link Vi, let us introduce 
(wo(j),wi(i), uj 2 (i) , u)3(i)} the corresponding base in a /owr-dimensional space, and 
consider </> a,,3 (i), a, (3 £ {0, 1, 2, 3}, the set of operators defined on this space by 



3.2 Chaos game representation 
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In particular (j) a,a {i) is the projection onto the state oj a at site i, the remainder 
of the sequence being left unmodified. With this notation, the generator takes the 
form 

g(i) = X+ [4> a+1 ' a (i)(f) a < a+1 (i + 1) - c f> a ' a (i)(f, a+1 ' a+1 (i + 1)] 

+ a- [ct> a ' a+1 (i)(i> a+1 > a (i + 1) - 4, a+1 ' a+1 (i)4, a ' a (i + 1)] 
+ 7 + [4> a+1 ' a (i)4> a - 1 > a (i + 1) - 4> a ' a (i)4> a+2 ' a+2 {i + 1)] 

+ 7 " [ ( /) a ' a +\i) ( t) a - a ~ 1 (i + l)-(f > a > a (i)(f > a + 2 > a + 2 (i + l)] , (3.1) 
and commutation rules write 

[4>^(i), <t^(j)} = % (<W a,!/ (*) - 5 av ^(i)) . 

The form ([Tl]) could be used to construct a transfer matrix representation 
For the moment, one can simply remark that the Lie algebra £ generated by the 
4>^ v is a Z2-graded Lie algebra [Q, so that it can be decomposed onto two Lie 
sub-algebras denoted by Cq and Ci, with the commutation relations 

{[Ab£o] C Co, 
[A, A] c Co. 

This leads to introduce the following linear combinations 

k p,i = (j> a+q < a , (3.2) 

ae{0, 1,2,3} 

for p,q € {0, 1, 2, 3}, with the commutation rules 

where the sums p + p' and g + g' are taken modulo 4. One easily verifies that, when 
q varies, K°' q and n 2 - q generate Co, while k 1 ' 9 and K 3q generate C\. In addition, 
from the very definition ([s"^), the quantities 



Qo 



















do commute with the generator. Later, by means of a more intuitive representation, 
we will be able to interpret the above formalism in terms of conservation of particles. 

3.2 Chaos game representation 

A way commonly used in DNA modeling (pH, Q and references therein) in order 
to visualize dynamics in the configuration space consists in encoding the random 
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3 Sequence coding and generators: two formal approaches 



walk as an array of integers at € {0,1,2,3}, i € N, and to associate with each 
sequence a complex number z, where 



N-l 



= ^ ex p(" n V 2 )- 



p=0 



Graphically, this leads to the successive inclusion of disjoint square sectors repre- 
senting the filtration defined by the random walk (see figure 3^ ) . To each random 
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Fig. 3.8: Chaos Game representation of the first two links of the chain (left) 
and change of scale (right) representing the next link. Arrows indi- 
cate the transitions 7 (center) and A (on the edges). 

walk of length N corresponds a square of area A~ N centered at the related point z. 
There are exactly 4^ elementary squares and it becomes possible to define a discrete 
probability measure n(z, t), constant on each such square, and whose evolution at 
time t is given by the forward Kolmogorov's equations 

~^(*'*) = / dyn(y,t)g(y,z). 

where g(z, z') denotes the kernel of the infinitesimal generator in the above rep- 
resentation. Let go be the part of the operator which operates at the upper level 
of the representation, that is on the first two links (see Figure |3~^) , and let D the 
dilation operator mapping z onto 2z. Then g writes 

g = g + DgoD- 1 + D 2 goD~ 2 + ... 

Moreover, as N — > 00, g must satisfy the consistency relation 

g = go + DgD- 1 , 

which can be split into two parts 

def 

9 = 9i +92, 

where 

' g\ = g + D 2 9l D- 2 
g2 = DgoD- 1 + D 2 g 2 D~ 2 
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It appears that 31,32 can be easily diagonalized, since 30 is diagonalizable and 
commutes with D 2 goD~ 2 and all higher order terms. Hence 31,32 can be endowed 
with a simple tensorial structure which, in matrix form, reads 

00 01 02 03 11 12 13 10 22 23 20 21 33 30 31 32 



^0 = 



0] 

[-A+] 
[-7] 
[-> 


[A 
[7-] 

-] 


1 


[7 + ] 
[A+] 


[7+] 
[A+] 


0] 

[-A+] 
[-7] 
[-> 


[A 

-] 


1 




[7+] 

[A+] 


0] 

[-A+] 
[-7] 
[-> 


[A 
[7-] 

~] 


[A 

[7-] 


1 


[7+] 
[A+] 


[0] 

[-A+] 
[-7] 
[-> 



The self similarity is revealed by the fact that each 4x4 block element becomes 
16 x 16 when one takes an additional link, in such a way that Pq is added to itself 
in these blocks. The problem is intricate since 31 and 32 do not commute. When 
7 + = 7~, Pq can be diagonalized and its eigenvalues are (9 times degenerated), 
— (A + +A~) (4 times degenerated), —7 (twice degenerated) and —2j (simple). These 
simple tensorial structure simply reflects the scale invariance in the chaos game 
representation or the translation invariance in the initial representation. When 
N — > 00, it should be possible to exploit this symmetry to build an iterative map, 
aiming at generating the invariant measure or the self-organized critical one in the 
glassy phase. Before doing this, we will propose another formulation of the problem 
in terms of particle hops with exclusion. 



4 Coupled exclusion models and thermodynamic limit 



4.1 The mapping 



In the last section, we used a representation helping to visualize the configura- 
tion space. There transitions between configurations were expressed as exchanges 
between square domains. We may go a bit farther on, using the fact that muta- 
tions are always either vertical or horizontal, and with opposite directions between 
both modified links (see Figure |4.9|) . This suggests indeed to recode each link j, 
j = 1, . . . , N by means of two binary components € {0, 1} and S {0, 1} thus 
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4 Coupled exclusion models and thermodynamic limit 



establishing a bijection aj — > (s",s*) such that 

0^ (0,0), 

1- (1,0), 

2 - (1,1), 

3- (0,1). 

In this scheme a transition on a link does touch only either of its components. 
Indeed the infinitesimal generator is the sum of two terms: the first one acts on 
the sequence {sf }, with rates conditioned by the sequence {s b } and vice-versa. We 
have thus a Markov process with state space 

{S a ,S b ) = {{s a 1 ,s b ),...,{s a N ,s» N )), 

a 2iV-dimensional boolean vector. Then, for an arbitrary function / : (S a , S ) — > C, 
the generator decomposes into 

N 

G = J2K(i) + h b (i), (4.1) 

i=l 

where 

M*)/(S Q , S b ) d = f A+(i)*?5? +1 [/ ((*?, Sl ), (0, a J), (1, s b +1 ), (s^r, 4)) ~ /(S a , S b )] 
+ (t)5? S ? +1 [/ ((*?, S J), (1, 4), (0, -, K, 4)) - /(s a , S h )] , 

h h (€)f(S\ S b ) d = f A+(i)*J*J +1 [/ {{sl, s\), «, 0), 1), {s%, s» N )) - f(S a , S b )] 
+ KHHsli [f ((*?, s b ), {s?, 1), {s? +1 , 0), k, 4)) - /(s a , s 6 )] , 

using the boolean notation s = 1 — s. In order to express these generators in terms 
of local operators, let us define the following pairs of operators (aj,Oj) and (bj,bj): 
they leave unchanged the sample path, but the first [resp. the second] component 
of link j , and they are given by 

a 3 f(S a ,S b ) d ^ f s?f{{s a 1 ,s b ),...,{0,s b ),...,{s a N ,s b N )), 

b J f(s a ,s b ) d = f 4/(K, s ?),...,(^,o),...,(^,4)), 

a] f{S a ,S b ) d ^ f ^f{{s a 1 ,s b ),...,{s^l),...,{s a N ,s b N )), 

b] f(S*,S b ) d ^ f s b f{{s a 1 ,s b ),...,{l,s b ),...,{s a N ,s b N )). (4.2) 

In terms of these operators, we also have (for (a) species) 

K{i) = Aj(i)(l - a ia \ +1 )a\a l+1 + A+(i)(l - a i+1 a\)a\ +l ai. 

Transition rates can be computed by inspecting the different cases, which yields 
the following expressions. 

'AJ-(i) = s b s b +1 \^ + 44 +1 A± + s b s b +il^ + 44+i7 ± , 

(4.3) 

Af (») = ^s» +1 A± + « +1 AT + s? s ? +l7 ± + s?s? +1 ^. 



4.1 The mapping 
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Fig. 4.9: Horizontal and vertical exchanges corresponding to transitions 
when link i has the value 0, according to the coding of section 



3.2. Rules for other transitions follow by rotational symmetry. 



The generator (4.1) represents two coupled systems of particles moving on a one- 
dimensional lattice with exclusion (i.e. there is at most one particle of each species 
per site fig)), 

These particles perform random elementary jumps to the left or to the right. Ob- 
viously, both systems are interlaced: the jump rates Aj (i) of species (a) at site i 
are conditioned by the stat es of particles of species b at site + and conversely 
according to relations (O). Setting 



A 



7 



A+ + A" 



and fi 



A^ 



A" 



def 7 



■7" 



and 5 



def 7 



' 7" 



2 2 
(with similar definitions for A a , Afc) the sub- generator of species (a) rewrites 

h a {i) = X a (i)(a\ - a\ +1 )(ai + i - ai) 



(4.4) 



+fj, a (i) a\{a i+ i - ai) + (a\ - a\ +1 )t 



where X a (i) and fi a (i) denote respectively the diffusion coefficient and the drift 
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4 Coupled exclusion models and thermodynamic limit 
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Fig. 4.10: Correspondence between chaos game representation of the random 
walk (a) and one-dimensional models of particle diffusion with ex- 
clusion. In (b), the sequence {s a } [resp. {s b }] determines the 
profile of the diffusion [resp. the distribution of particles], drawn 
for 7 1 * 1 = X^; in figure (c) the role of the particles has been ex- 
changed. 



(with the same notation as in (4.4)). After a routine algebra, we obtain 



A Q « = A + ( 7 -A)( Ss b +1 - Sl b ) , 

H a {i) = M (l-4-^ +1 ) + «5(4 +1 -4), 

A 6 « = A + ( 7 -A)« +1 -<) 2 , 

M&« = K S ? + < +1 -l)+5«-< +1 ). (4.5) 

It is worth remarking A a and A& are strictly positive except when 7 equals zero, in 
which case they might be zero at certain points. The interaction term takes the 
form 

V a (i) = 2\ a (i)s1s1 +l + 8t(jia{i) - Va{l - 1)), 

where ^ a {i) can be interpreted as a potential. A convenient representation of the 
system is to draw a one-dimensional profile from the sequence {s b } (positive or 
negative slope depending on whether s b equals or 1), t o sketc h the probabilistic 



inclination to turn left or right, as shown in Figures |4.10| b, |4.1C| c. 
In the particular case 7 ± = A* , we get 

'A Q ± « = i(A++A-)±(,s b -i)(A+-A-), 

(4.6) 

A±(i) = i(A+ + A-)T«-|)(A+-A-), 



4.1 The mapping 
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and the sub-generator of species (a) takes the form 

h a (i) = Xa\(a l+1 + 04-1 - 2a,) + ll{2s\ - l)(af - aj+i) (a* + a i+1 ) + V a (i), 
where V a (i) is a diagonal term, 

In this particular case, A represents explicitly the diffusion constant of an isolated 
particle, with drift fi, the sign of which is determined by s\. V a (i) is reminiscent of 
the interaction between particles coming from the exclusion constraint. 





Fig. 4.11: (a) Elementary transition due to the jump of a type (a) particle, 
and the corresponding deformation of the profile defined by s a (b) 

( 7 ± = A±) 

The distribution of particles labeled (a), submitted to the diffusion defined this 
way, is given by the sequence {s a } (0 or 1 particle depending on the individual site 
values of s Q ). A complementary model is obtained by exchanging the roles of s a and 
s b . In addition, elementary transitions of the system correspond to jumps in the 
left or right direction of particles (a) and (b). In the complementary formulation, 
these transitions are materialized through modifications of the profile determined 



by s a or s b (see figure |4.1l| ). Therefore, from this viewpoint, we can formulate 
the dynamic of subsystem (b) in terms of a KPZ model |l7j |, in which the noise is 
produced by the distribution of particles (a). With this formulation, the conserved 
quantities (pointed out earlier) can be obtained in a straightforward manner, since 
they simply express conservation of particles. If boundary conditions are such that 
particles cannot escape from the system, the population of both species is conserved. 
This is for example the case when we impose periodic boundary conditions or also 
when extremities of the chain are fixed. 

Suppose we fix the total amount n a and rib of particles (a) and (6), this results then 
on the random walk by the fact that no + n^ = n a and no + n\ = rib are fixed (rii is 
the number of links i) . Since no + n\ + ri2 + m = N it then easy to convince oneself 
that this is equivalent to fix no — ri2 and m — ri3, which is enough to determine the 
respective positions of the initial and final points of the walk. 

In addition, except when 7 = 0, we get an easy way to determine the irreducible 
classes of the system. In fact, for closed systems with a fixed number of parti- 
cles, once the population of each species is fixed, just start from the configuration 
(11. .1100. .00) for both species where all particles have been disposed to the left. 
Then, owing to possible consecutive jumps to the right, which for 7 7^ are al- 
ways authorized, one can reach any arbitrary configuration. This shows that for 
closed systems irreducible classes are indexed by the number of particles in each 
species, which corresponds to the separation between extremities of the walk. On 
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2 




(c) (d) 

Fig. 4.12: (a) Stable configuration for closed chains, (b) Corresponding rep- 
resentation in terms of exclusion models, (c) Glassy state with 
N = 5000. (d) Corresponding KPZ landscape with the density of 
trapped particles represented in white. 



the other hand, for open sample paths with free boundary conditions, irreducibility 
holds as long as particles can both enter and leave the system. An interesting point 
is the way stable configurations, which are numerically observed, are represented 
by me ans o f this exclusion process formulation. In agreement with the intuition, 



Figure |4,12| depicts a stable situation where particles are trapped in a well. As 
this re main s true with the complementary representation (see the lower part of 
Figure 4. 12j b) , the following iterative scheme could be used to generate the global 



invariant measure: let particles (a) evolve assuming the dynamics of particles (b) is 
frozen; then, once this conditional stationary regime is reached, switch to particles 
(b) conditioned by particles (a), etc. Translated mathematically, we proposed the 
following iterative scheme in order to capture the invariant measure 

X n+1 = lim E (S a {t) | S b (t) = Y n , S a {0) = X n ) , 

t-^-oo ' 

Y n+1 = lim E (S\t) | S a {t) = X n , S b {0) = Y n ) . 



4.2 Conditional equilibrium 
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Numerically, the sequences of random variables X n and Y n seemingly converge for 
77 < 50 to the stationary variables S a (oo) and S b (oo). 

Up to an abuse of notation, we shall often identify the random process with any of 
its sample paths. For instance we shall simply write P(S a \ S b ) [resp. P(S b | S a )] 
the conditional invariant measure of particles (a) [resp. (b)] when the dynamics of 
particles (b) is frozen. Then the iterative scheme can be reformulated as 

P n+1 (S a ) = p ( sa I S b )Q n {S b ) 

{S>>} 

Q n+ i(S b ) = ^2 P{S b I S a )P n (S a ) 

{S°} 

where P n [resp. Q n ] represents the probability measure of the S a [resp. S b ] after 
n steps. If this iteration converges, the joint probability for the invariant measure 
will be given by 

P(S a , S b ) = P(S a I S b )Q(S b ) = P(S b I S a )P(S a ) 
in case these two expressions coincide. 



4.2 Conditional equilibrium 
4.2.1 The case of a stretched walk 



Let us have a look to a special case which can be solved exactly. It will provide 
some hints about scaling of the parameters when N — > 00. 

Consider a random walk with fixed extremities, and consisting only of links oriented 
either to the north or to the east, i.e (at € {0, 1}, i = 1 . . . N). Here folds do not 




Fig. 4.13: Stretched random walk and asymmetric exclusion process. 

exist (see figure ( 4.13| )), so that s b = everywhere and solely transitions A ± take 
place. It turns out that the invariant measure has the product-form 

P {Oil, ■ ■ • , Q-n) = P ai Pa 2 ---Pa N - 
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Indeed, introducing the occupation rate of particles (a) qi = pi(i) — 1 — Po(i), we 
have the balance equations 

(4.7) 



Setting 

we obtain a geometric series 



A + (l - qi)qi+i = A <&(1 - q i+ i). 

def Qi 



1-Qi 



In 



the ratio 





I 






)- 


r exp 


V 


A 


f - A" 


N 


A 


f + A- 



ilog 



1 + 



(4.8) 



A' 



already introduced in (2.1) and (4.4), gives the typical scale n c ~ above which we 
obtain straight aligned patterns. Here ro is a normalizing constant, which permits 
to fix the expected value Dq of the particle density 



that is 



1 N 

i=l 



»=1 " j=l 

Letting iV — > 00, we analyze the limiting process under the scaling, 

N^oa A ' 

where 77, up to some abuse in the notation, is now a parameter independent of N. 
This is tantamount to assume that fx is implicitly a function of N. 



def ,. 

77 = lim 



(4.9) 



Fixing x = i/N and taking the expansion with respect to r]/N in (|4.q), we get the 
limit equation 

r(x) = r exp(— 2rjx), (4-10) 

which implies in turn 



D a = I d.r j 1 - 
Consequently, 



1 



1 + r cxp(— 2-qx) 

smh{rjD ) 



— lo 1 + r ° 

277 1 + r cxp(-2?7) ' 



sinh?7(l — Dq) 



The integral of the particle density, taken as a function of a;, is shown in Figure [4.13 
and is given by 



h(x) = du 1 - 



1 



1 



log 



l + r 



1 + ro cxp(— 2rya;) 



1 + ro exp(— 2?7u) 

This asymmetric exclusion model can be solved under more general conditions, in 
particular with open boundaries, using matrix methods plL The above simple 
example confirms the observed fact that the correct scaling parameter is indeed 77, 
and also somehow explains why the chain remains Brownian when 77 < 1. 



4.3 Closed curves: weak convergence for large N 
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4.2.2 The general case 

Once the sequence {s b } is given, particles (a) form a unidimensional exclusion 
process in an inhomogeneous environment. Its transition rates to the right or to 



the left at position i have been given in (4.3). 



Our basic claim is that the reasoning of section 4.2.1 still holds for a closed system 
where the number of particles is kept constant. Exactly this means that, as long 
as there is no current in the system, conditional detailed balance equations of type 



(4.7) are still valid at steady state, just replacing A* by A„ (i). 



We present no proof of this fact (which could likely be obtained by coupling argu- 
ments), and write brutally the equilibrium equations 

A+(i)(l - «?)«?+! = A"(*k a (l - q« +1 ), i = l...N, 

where sites N + 1 and 1 are identified for periodic boundary conditions, and qf 
[resp. q\\ is the random variable equal to the conditional probability of having one 
particle of type (a) [resp. (b)] in position i given the sequence {s b } [resp. {s a }], 
that is 

q? = E[sl\S b ], 9 ? = E[ S »|S»]. (4.11) 

Setting 



l-qt 
we obtain 

log[r? +1 ] - Iog[rf] = lo g(^||) • ( 4 - 12 ) 
By using ( |4.3| ), an easy algebra based on the boolean character of the s\ yields 

log 5$ = (1 - 2s b ) log ^ + (s?+1 _ s >) log ^T^ + fj , (4.13) 
XJ(i) A + /i (7 + d)(A-/i) 

with a similar equation for type b particles. In addition we observe that the con- 
straints 

r a _ a b _ b 

r N+l — 'D ' N+l — "l) 

for periodic boundary conditions, will be automatically fulfilled as long as the sys- 
tem of particles (b) is globally neutral (i.e. particles and holes have the same 
cardinality), in which case 

N 



4.3 Closed curves: weak convergence for large TV 

We will now combine the stationary product forms obtained for each particle 



species, according to the iterative scheme proposed at the end of section |4.1 
Thr ough out this section, the dependence on N of the random variables q^ : q\, given 
by ( |l.ll| ), is kept implicit for the sake of shortness in the notation. 
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4 Coupled exclusion models and thermodynamic limit 



Introduce {wf,i > 0} and {w b ,i > 0}, two families of independent and identically 
distributed Bernoulli random variables with parameter 1/2 and values {1,-1}. 
Without further comment, we assume that (S a , S b ), {wf , i > 0},{w b ,i > 0} are 
defined on the same probability space. 

Lemma 4.1. Let au,k > 1, denote a sequence of complex numbers satisfying the 
condition sup fe \otk\ < oo. 

There exists a probability space such that 

N JV 



Jf E a * 8 l = E + °Wk) + 0(N- 2 ) a.s. 

fc=i fe=i 

]V E a ^ = n E + + 0(N- 2 ) a.s.. 



fc=i 



fe=i 



where 



(4.14) 



(4.15) 



Proof. The starting point is the straightforward decomposition of the probability 
measure P(S a ,S b ) as 

P(S a ,S b ) = P(S a I S b )P(S b ) = P(S b I S a )P(S a ). 

Considering the invariant measure of the process (S a (t), S b (t)), we claim the condi- 
tional probabilities P(S a | S b ) and P(S b | S a ) coincide with the equilibrium proba- 
bilities obtained in the last section for each particle species. In the present situation, 
this is tantamount to writing the equation 



JY 



P(S a \S b )=Y[(s?q? + s?(l- q ?)), 



where the qf 's depend implicitly of S b . To analyze more precisely the coupling 
between the two families, we introduce the Laplace transforms 



J 1 ^ * 



k=l 



N 



n[ i+ ^( e ' 

k=l 



Then 



exp 



N N 

N E °M2 + 2N* E - ) + °(^ } 



fc=i fc=i 
On the other hand, starting from the equality 



exp 



1 ™ 1 1 r c 1 w N a 

- <**(«? + «) = E exp - ]T a fc8 g + £ log cosh Z&L 

k=l ' -I L fe=l fe=l 



with regard to the (a) species, we observe that the value of er£ given in (4.15) gives 
the matching 



( i N ii r r i w 

exp 1 }v E afcS fc f = Eexpj — 2^«fc(«fc + < 7 />fc 
fc=i J J L I fc=1 



(4.16) 



4.3 Closed curves: weak convergence for large N 
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Since all rand om variables at stake are uniformly bounded, equation ( 4.16; ) yields 
at once ( 4.14 ), but only in distribution. To conclude the proof of the lemma, we 
make use of transfer and coupling theorems due to Skorohod and Dudley (see , 
theorem 4.30 and corollaries 6.11, 6.12), which allow to switch from equalities in 
distribution to almost sure properties, since on the original probability space the 
right member of system (4.14) is a measurable mapping of the left one. ■ 



4.3.1 Fundamental scaling, thermodynamic limit and fluctuations 

For any i, 1 < i < N, we put ad libitum x — i/N, < x < 1. 
Proposition 4.2. Under the fundamental scaling 



' A* V , ( 1 A 



the weak limits 



6 - = o(L 

7 \N 



p a (x) = lim q a xN and p b (x) = lim q b x 



N 



N- 



N 



(4.17) 



exist and satisfy the autonomous system of deterministic nonlinear differential equa- 
tions 

dp a (x) 



dx 



A V p a (x)a~p a (x))(p b (x)-l), 



(4.18) 



^M = -4 V p b (x)(l-p b (x))(p a (x)- 1 -] 
In addition, the assumed closure of the original random walk imposes the relations 

Jo P a {x)dx = J p b (x)dx = |, 

p a (x + l) =p a (x), (4-19) 
^p b (x + l) = p b (x). 



Proof. Taking the expansion with respect to N in equations (4.12) and (4.13), and 
using lemma iA , we get after some algebra 



fc-i 



1 i=i 



(4.20) 



Since N is a parameter and x rather stands for a variable, it is be convenient to 
introduce the following functions of x 

=^ xN] , r%(x) d = f ^ *%(x) d = f J q « N (x)(l-q%(x)), 

1 Hn\ X ) 
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and similarly for the (b) species. 



To omit some tedious technicalities, we will only sketch the remaining lines of the 
proof. 



First, it is not difficult to see that by restricting the expansion in (|4.16|) up to terms 
of order iV -1 , and using the definition of r£,r£, we come to the simplified system 



k-l 



-1) 



3 = 1 



k-l 



+ 



1 + r? exp [^£(2^-1) 



r\ exp 



4 = 



3 = 1 
3 = 1 



k-l 



3 = 1 



a.s., 



(4.21) 



In a second step, it can be shown from ( 4.21 ), as in a purely deterministic context, 
that the quantities q^ix), q N (x), form Cauchy seque nces , hence converging, f or al l 
< x < 1. To see that the deterministic limits ( 4.17 ) exist and satisfy ( 4.18 ) 
is straightforward by approximating discrete sums by Riemann's integrals. This 
yields the differential system 

which has similarities with the famous Lotka-Volterra equations, where x plays here 
the role of the time. It is worth noting that ( |4. 1C| ) is immediately revisited, taking 
merely p b (x) = (i.e. the density of particles (b) is kept constant). ■ 



It is also tempting to get some insight into fluctuations around the above deter- 
ministic limit. This will be achieved by establishing the forthcoming central limit 
theorem. 



Proposition 4.3. Under the fundamental scaling, the weak limits 



g a {x) = lim VN 



q%{x)-p a {x) 



g°(x) = lim VN 



q b N (x)-p b (x) 



exist and satisfy the system of stochastic differential equations 
dg a {x)=Arj\(j 2 a {x)g b {x)dx + a b (x)dW b (x) 



where 



dg b (x) = -4?7 al(x)g a {x)dx + a a (x)dW a {x) 



a a (x)^ ^/p-(x)(l-p-(x)), a b (x) d ^^pb(x)(l-pb(x)). 



(4.22) 



(4.23) 



4.3 Closed curves: weak convergence for large N 
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Proof. Letting 



[xN] [xN] 

X« N (x) ^ -j= £ w%, X b N (x) d ^ f J= £ «*> 



we consider the two white noise processes 



W a (x) = lim X a N (x), W\x) = lim X b N (x). 



N- 



N- 



We call on a strong approximation theorem, which is a refinement of the invari- 
ance principle (see |18|: it shows how to construct W a (x) and W b (x) on the same 



probability space as X^ and X^) in such a way that 



sup \X N (x)-W(x)\=0( l ^£) 
o<z<i v VN ' 



Arguing as in the derivation of proposition 4.2 



we can write 



log 



r b N (x) 



du(2q b N (x) - 1) + -% f a b N {u)dW b (u) + o(^jf 



r b N (0) 



-2r,£ du(2q a N (x) - 1) - i| £ a%(u)dW a (u) + o(^) , 



(4.24) 

where the stochastic integrals are taken in the Ito sense (see J2^|). Defining g%{x) 
and g%(x) by the equations 



p a (x)+a 2 a (x)g%(x), 
P b (x)+a 2 b (x)g b N (x), 



then putting these expressions into (4.24) and differentiating (4.24) (details are 
omitted), we are lead to (4.24. The proof of proposition 4.3 is completed. ■ 



Remark Let us comment on the scaling of S. Cont rary to the scaling of p, which 
is naturally dictated by homogeneity (the sums in ( 4.20 ) remain meaningful after 
dividing by N, when N — > oo), it is a purely dynamical consideration which dic- 
tates the scaling of 6, Indeed, fi and 5 are associated to time constants r M and 
Ts- The quantity t m represents the typical unit of time for a free particle (a) or 
(b) to drift along the system over a finite distance, whereas ts is a time-scale for 
rotations of vertical or horizontal fold M2 of the chain, remembering that S stands 
for the detuning between 7 + and j~ defined in section H. Therefore either these 
time-scale are coherent and 5 is rescaled, otherwise rotational motions of M2 occur 
at a shorter time scale and a different analysis has to be conducted, since motifs M2 
(which correlate hole-particle pairs of species (a) and (b)) reach their equilibrium 
distribution before particles have enough time to move along the system. Numer- 
ically we could not perceive any specific effect related to S, so that we restricted 
ourselves to the above fundamental scaling. 
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4.3.2 Second order phase transition 



Here we focus on on the deterministic ( |4.18| ) part of the equations. Using the 
notation 

v a (a) d = 2p a (x) - 1 , v b (x) d = f 2p b (x)-l, 



apart from the trivial solution v a (x) = Vb(x) = 0, we obtain from fll.lq ) 

f a (x) 9^ Q (x) _ t%(a;) cH(:r) 



1 — f„(a;) <9x 1 — v%(x) dx 



or, after integration, 



< C < 1, 



since |i/ | < 1 and \v b \ < 1. Plugging the last relation into (LIS) leads to 



dv a (x) 
dx 



r, 2 [l-ul(x)][l-C-^(x)}, 



(4.25) 



the solution of which is the standard Jacobi elliptic function 

1 



v a (x) 



VT^C 



sn(i]x, vT— C). 



Finding the constant C is equivalent to compute the fundamental period of these 
functions. Hence, denoting by X(C) the period of v a {x), we have solve 



X(C) = 1. 



(4.26) 



From ( 4.25 ) we see that v a {x) is bounded by — \f\ — C and \J\ — C. Taking into 
account the constraint (J Q v a {x)dx = 0), X(C) is exactly given by 

x{C) = VT^c) = - t t 
v ' v v 2' 1 nJo ^[i_^][i_(i_c)^]' 

where F is the standard elliptic integral of first kind. 

X{C) is a decreasing function of C on ]0, 1], reaching its minimum for C = 1, so 
that 

X(C)>X(1) = -. 
Thus appears is a critical value for rj, namely 

r]c = 2n. 



When -q c < 2ir, ( 4.26 ) cannot be fulfilled and we are left with the trivial solution. 

When r] c > 2ir, it is straightforward to compute the arithmetic area Si(rj) of the 
first winding sector, since it is indeed the only non- vanishing sector. Setting 

1 f u 

h x {u) =2 dv[v a (v) + v b {v)]dv, 

1 f u 

h y( u ) = j / dv[v a (v) - v b (v)\dv, 



4.3 Closed curves: weak convergence for large N 
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Si(t]) is simply the area enclosed by the curve (h x (u),h y (u)) u € [0,1], which is 
given by 



du 



h dh v 

n x a 

ou 



dh x 
' du 



or, after some algebra, 



V /(1_^2 )(1 _ C _^2 ) 



log 



rl 



1 -v 



keeping in mind that C is also a function of rj. The corresponding curve displayed 
in figure 2^.b matches pretty nicely all numerical observations, in particular with 
regard to the critical value r] c . 



4.3.3 More about fluctuations for rj < r\ c 



When rj is under the threshold r) c , the deterministic part becomes trivial, and we 
are left with fluctuations. This corresponds basically to the observations shown in 
Figures 2J5.& and ^.b. Inserting 



in ( |4.23| ), and setting 

g{x) = g a (x) + ig b (x), W(x) = W a (x) + iW b (x), 
we get solutions of the form 

g(x) = -2irj f e" l(u - x) dW{u), 
Jo 

which corresponds to 



(4.27) 



<» = 5 + Jv 



sin 



(r)(u- x))dW a {u) + cos(r?(u- x))dW b (u)) + dW a (x) 



^ = l + 7N 



cos 



(ri(u-x))dW a (u)~sm(i 1 (u-x))dW b (u)) + dW b (x) 



Hence we have derived an equivalent process, which up to order 0(^7=), describes 
the curves numerically observed. Letting 



[xN] 



3=1 

[xN] 



the so-called equivalence says precisely 

dh%(x) =h%(x + ±)- h%(x) = 2q%{x) + -j=a%(x)dW a (x) - 1 + o(-j=) 



28 



4 Coupled exclusion models and thermodynamic limit 



and the same holds for dh b N (x). Introducing the complex function 



}in{x) = h a N (x) + ih h N (x), 



we have 



,m i \ 2ir ) 
atiN(x) = — 



which, after integrating by parts, yields 

h N {x) = -L [ / 2e l1l{u - x) dW{u) - W{x)\ + of-L) . 

V N k/o J ^ \ N ' 

At this point it is possible to reconstruct the curves observed numerically, remem- 
bering that the discrete displacements (dxi, dyi) in the plane are expressed in terms 
of sf and s° as 

U Xl = 1-8* -8%, 

1 dyi = sf- si 



According to (4.27), we define 

[xN] 



h(x) = lim —== (dxk + idyk) and Z(x) 

VJV fc=0 



i - 1 



W{x). 



Then we have 



h{x) 



2 exp[iry(u — x)]dZ(u) — Z[x). 



This above equation accounts for the windings of the Brownian curve observed in 
figure 2.3.b. When — > 0, h{x) coincides with the standard Brownian motion Z . 



4.4 Burgers equations in the fluid limit 

In this section with give a formal derivation of dynamical equations describing 



the fluid limit, which includes the steady state solutions obtained in |4.3| , without 
pretending to a rigourous presentation. 

We start from the heuristic assumption that the conditional independence of the 
sf [resp. sf], which is realized at time t = and at steady state, remains valid for 
all fixed time t, up to order 0(N~ 1 ). Exact proofs of this fact could be provided 
by adapting (up to sharp technicalities) some lines of argument proposed e.g. in 



1 19, 10]), together with a mean field type approach for the convergence of the semi- 



groups of the underlying Markov processes indexed by N. 

Considering the stochastic variable expressing the current of particles (a) between 
sites i and i + 1 at time t 



tf(f) = X+(i,t)sf(t)sf +1 (t)-X-(i,t)sf +1 (t)sf(t), 



<p\{t) = A+(M)«?W4fi(*)-A 6 _ (M)4fi(*)S<(*)> 



4.4 Burgers equations in the fluid limit 
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we define the conditional expectation 

J? (t) = E[tf(t) | S b ], J$(t) = | S a ]. 

On account of the particle conservation principle, we have locally 

^E[<(t)] + E[^(t)-^_i(*)] =0- 



(4.28) 



Then introducing the time dependent expectations q k (t) as functions of the sample 
path S b 

q a k (t) = E[s a k (t) | S% q b k (t) = E[s b (t) | S a ], (4.29) 



we can write by ( 4.28 



-Y,a k (q a k (t) - q a k (0) + / dr(j a k (t) - J^(t))) = 0, (4.30) 

k=l •'° 

with again the condition sup fc |ctfc| < oo. Expressing the almost conditional inde- 
pendence of the s k 's, we have 



A' 



- ^a fc E t A+(M)tf(*)(l - qt +l (t)) - X~(k,t)q k+1 (t)(l - q a k (t)) 



k=l 



1 N 1 

_X;a fc JJ(*) + o(-). 



k=l 



To be consistent with the procedure developed for the stationary regime case, we 
substitute to the S b the so-called equivalent set {q\ + cr k w k } into the expressions 
of the rates given by 



logAj(i) 



logA + 2(^ + 4 fl -2 s ^ + i) 1 c.g 



2 

A 



±l(l-s b -s b +1 )± 6 -(s b -s b + 



according to the fundamental scaling. This insures that A+ and A~ remain perfectly 
correlated. As for the the fluid limit, the procedure a moun ts simply to replace s b by 
q b in A„ (i), and to approximate all discrete sums in ( 4.30 ) by Riemann's integrals, 
for arbitrary a(x). This yields the continuity equation 



dp a (x,t) _ dJ a (x,t) 



dt 



dx 



(4.31) 



where p a {x,t) and J a {x,t) are the deterministic continuous counterparts of q k (t) 
and J k (t). The current is now given by 



J a (x,t) = D 



1 atl a\ft o b\ dp a (x,t) 

2r]p (1 - p )(1 - 2p J — — 



exp(2p 6 (l- /9 6 )log^J, (4.32) 



where we have introduced the diffusion constant 
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5 Conclusion 



A being implicitly taken as a function of N. This scaling is confirmed numerically 
(see figure ||l|c). Actually we observe that the parameter 7 controls the dynamics 
of the system through the definition of an effective diffusion constant. For particles 
(a), we have 



D a (x,t) = Dexp 2p b (x,t)(l- p\x,t))\oi 



with the corresponding relation for particles (b)). In the particular case 7^0, 
this constant vanishes, except at loci where the density of particles (b) has no 
fluctuations, that is a 2 — p b (l — p b ) = 0), in which case D a {x,t) = D, as to 
be expected from the analysis of the stretched walk. When 7 = A, we obtain a 
dynamical system of deterministic equations 



,„ - ^^-2D,|;[p.(W..>(l-2/)]<z, t ), 



These equations belong to Burgers' class. When taking one of the two density 
species (say (a)) to be a constant p a = or p a = 1, the density of particles (b) 
is driven by an ordinary Burgers equation, describing the evolution of a stretched 
walk. For an arbitrary 7, the steady state solution of (4.31) is tantamount to let the 
current vanish in (4.32), which after integration gives system ( 4.1g| ) independent of 
7, as to be expected. 



5 Conclusion 



The model of the discrete event system presented in this report turned out very 
friendly for simulation runs. Although the dimension be small (curves in the 
Euclidean plane), several basic phase-transition phenomena have been observed, 
among which a glassy phase. As we stove to point out, there are many ways to 
describe this system, which bring to light connections between various stochastic 
and algebraic formalisms. 

Nonetheless, in our opinion, the most efficient way toward concrete mathematical 
and physical properties appears to be on the track of coupled exclusion processes. 
This mapping allows also to address the continuous limit considered in the last 
section, and this approach can be a method for coding numerical simulations. In 
this manner, we have been able to observe that when we alternatively freeze one 
of the subsystems, letting the other one reach its equilibrium, the whole process 
attains its stationary state, and moreover much quicker. 

With regard to the dynamic, deeper investigations are needed in order to include 
fluctuations directly into the equations, and to clarify the role of the parameter 
S. This would give a firm starting point to get an insight into slow dynamics and 
into the non-linear excitations which are swarming in the glassy phase (metastable 
states). Open boundary conditions and presence of currents would also be worth 
investigating. Actually, it could be most rewarding (and this not hopeless) to 
generalize the model to higher dimensions and to find related concrete applications, 
for example in biology (evolution of RN A and proteins) . 
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